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ABSTRACT 

We investigate the cosmic evolution of low luminosity (ii.4GHz < 
lO^^W Hz^^sr"^) radio sources in the XMM Large Scale Structure survey field (XMM- 
LSS). We match low frequency selected (610 MHz) radio sources in the XMM-LSS field 
with near infrared /C-band observations over the same field from the UKIRT Infrared 
Deep Sky Survey. We use both the mean V/Vma.x statistic and the radio luminosity 
function of these matched sources to quantify the evolution of the co-moving space 
density of the low luminosity radio sources in our sample. Our results indicate that 
the low luminosity sources evolve differently to their high luminosity counterparts 
out to a redshift of z^O.8. The derived luminosity function is consistent with an in- 
crease in the co- moving space density of low luminosity sources by a factor of ^1.5 at 
z=0.8. We show that the use of the K — z diagram for the radio source population, 
although coarser than a full photometric redshift analysis, produces consistent results 
with previous studies using ~> 10 band photometry. This offers a promising method 
for conducting similar analyses over the whole sky with future near- and mid-infrared 
surveys. 
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1 INTRODUCTION 

The strong cosmic evolution of the most powerful radio 
sources was first deduced by iLongaij l|l966l ) from low fre- 
quency radio source counts. This evolution has since been 
confirmed by a number of investigations which imply that 
the co-moving space density of high luminosity radio sources 
has decreased by a factor of ap proximat ely 1000 between 
z ^ 2 and z ^ (e.g. iLain g. Rilc v fc Longairl 1 19831 : 
iDunlop fc Peacock|[l990i : IWillott et al.ll200ll). Beyond z ^ 

2 the evolution remains uncertain (e.g. Jar vis fc Rawlingg 
I2OO0I : IJarvis et all 



2001bl) bu t app ears to undergo a grad- 



ual decline (e, 
low-power radio sources is 
early studies of the radio 
the low luminosity 



WaU et all l2005h . The evolution of the 
less well constrained, but 
source counts indicated that 
could not be evol ving 



strongly as their high luminosity counterpa rts (jLongairl 
ll966l : lDoroshkevich. Longah fc Zeldovichll970l b Many mod- 
els of the evolution of radio sources thus divide the radio 
source population into two independently evolving compo- 
nents, a strongly evolving high luminosity component and 
a lo wer luminosity component with much wea ker evolution 
(e.g. I Jackson fc miill 19991 : IWillott et ablbOOlh . 



Radio galaxies can be classified into two mor- 
phologically distinct groups identified as Fanaroff- Riley 
class I (FRI) and Fanaroff-Riley class II (FRIT) sources 
(Fa naroff fc RilevllT974 ) . The FRI sources have higher sur- 
face brightness close to the centre of their radio lobes 
whereas the FRII sources have highly coUimated large-scale 
jets and bright emission hotspots at the edges of their radio 
lobes. The FRII sources are typically more luminous than 
the FRI sources with the division in luminosity falling at 
roughly I/i.4GHz = lO^'^ W Hz"^sr"\ This dividing luminos- 
ity may be dependent on the optical luminosity of the host 
galaxy, a higher optical host l uminosity results in a higher 
FRI/FRII division luminosity (|Ledlow fc OwenI 19961 ). It has 
been suggested that the FRI and FRII sources might respec- 
tively be associated with the slowly and rapidly evolving 
components of t he radio sour ce population (jjackson fc Walll 
IT999). However iRigbv et al.l (2007) find evidence that high 
luminosity FRI sources evolve as rapi dly as FRII sources of 
compa rable radio power. Furthermore, iGendre. Best fc Walll 
(|201Cf ) performed a detailed comparison of the radio lumi- 
nosity function of FRI and FRII sources at a number of 
redshifts. This analysis revealed that both populations expe- 



2 Mc Alpine & Jarvis 



rience luminosity dependent number density enhancements 
at higher redshifts, z~0.8~2.5, and that there were no sig- 
nificant differences between the enhancements measured for 
these two populations. These results would seem to suggest 
that both types of FR sources experience a common evo- 
lutionary history and thus cannot fully account for the ob- 
served dichotomy in cosmic evolution. 

Recent evidence suggests that a more important di- 
vision in the radio source population, distinct from the 
morphological classification, may be related to different 
modes of black hole accretion. Radio galaxies have also 
been classified according to the presence or absence of nar- 
row high-excitation em ission lines in the spe c tra of their 
optical host galaxies pfine fc Longairl 19791: iLaing et al.l 
11994 iJackson fc Rawlind Il997l : IWillott e^al]l2001^ . Ob^ 
jects without high-excitation emission lines are referred to 
as low-excitation radio galaxies (LERG). Low luminosity 
FRI galaxies are predominantly LERG's and most objects 
with high excitation emission lines are associated with the 
more powerful FRII sources. However the relationship be- 
tween FR class and the emission line classification scheme 
is not one-to-one as many FRII galaxi es have been foun d 
to be low-excitation radio galaxies (e.g. lEvans et al]|2006l ). 
It has been argued that these two groups correspond to 
different AGN phenomena powered by fundamentally dif- 
ferent modes of accretion. The LERG sources are powered 
by a radiatively inefficient accretion of the hot gas in the 
intergalactic medium, referred to as "radio" mode accre- 
tion. While the HERG are powered by radiatively efficient 
accre tion of cold gas, referred to as quasar mode accre- 
tion jEvans et al. I2OO6: Ha rdcastle. Evans fc Crostonll2006l . 
I2OO7I : [Herbert et al.. ,2010, ). iSmolcic et all (|2009l ') argue that 
the observed bimodal evolution of radio sources may be 
caused by differences in the evolution of these black hole 
accretion modes. 

Understanding the relationships between different 
classes of radio galaxies, and their cosmic evolution, is im- 
portant for understanding models of galaxy formation and 
evolution. Semi-analytic models of hierarchical galaxy for- 
mation predict that galaxies in the local universe should 
be more massive and have higher ra tes of star formation 
than is observed (| Bower et al.l I2OO6I ) . There is increasing 
observational evidence that radio activity from AGN may 
disrupt cooling flows in luminous early-type galaxies thus 
slowing or preventin g significant accretion of gas and fur- 
ther star formation IIFabian et al. I I2OO3I : iBirzan et af]|2004 
iRawlings fc Jarvij 120041 '). This negative AGN feedback has 
been successfully incorporated into models of galaxy forma- 
tion which are then able to better reproduce several fea- 
tures of the observed universe including the exponential 
cut-off in the bright end of the galaxy l uminosity function 
l|Croton et al.l I2OO6I : iBower et al]|2006l ). iBest et all (|2006l ') 
implied that the low luminosity radio sources may be the 
predominant contributors to this negative feedback effect. 

Current studies of the cosmic evolution of low power 
radio sources seem to imply that these sources are not 
evolving at all or experience only mild negative evolution. 
IClewlev fc Jarvi"3 (2004) found that the low luminosity radio 
sources with I/325MHZ < 10^^ W Hz~^sr~^ exhib i t no e volu- 
tion out to redshift z ~ 0.8, IWaddington et all (|200ll ') also 
found no evidence of evolution in their study of 72 sources. 
ISadler et al.l (|2007l ) find evidence of mild evolution out to 



z ~ 0.7 in their comparison of the radio luminosity function 
of sources in the 2dF-SDSS Luminous Red Galaxy (2SLAQ) 
and QSO survey with that of sources in the 6dF Galaxy sur- 
vey (6dFGS). This was consistent with pure luminosity evo- 
lution of the form (1 -f which ruled out the no evo- 
lution scenario at the 6a level. iDonoso. Best fc KauffmannI 
(2009) also find evidence of mild positive evolution in the 
z = 0.1 ~ 0.55 redshif t range. Using radio sources in the 
VLA-COSMOS survey ISmolcic et alT (,2009) find evidence 
that this mild evolution continues out to z ~ 1.3. All of these 
studies imply that the level of evolution in the low power ra- 
dio sources is signiflcantly less than that taking place in the 
high luminosity sources. 

In this study we use low frequency radio sources de- 
tected in the XMM-Large Scale Structure survey fleld to 
investigate the co-moving space density of low-luminosity 
radio sources. The low-frequency selection is preferred over 
high-frequency (e.g. > 1.4 GHz), as it provides an orienta- 
tion independent selection, as the low-frequency detects the 
optically-thin lobe emission, whereas high-frequency surveys 
contain a higher fraction of pol e-on sources where the opti- 
cally thick core dominates (e.g. Ijarvis fc McLure|[200^ ). 

In section 2 and 3 we discuss the radio and optical obser- 
vations used in our analysis. Section 4 outlines the method 
we used to cross match the radio sources with their optical 
counterparts. In section 5 we describe the method used to 
obtain redshifts for our sample. In the section 6 and 7 we 
discuss the y/Vmax statistic used to quantify the evolution 
of our sources. In section 8 we determine the radio luminos- 
ity function of the low luminosity sources in our observed 
field and use it to investigate their evolution. We present 
our conclusions in section 9. 

We assume the Hq—7Q km.s~'^.Mpc~^ and a f2A/=0.3 
and Q,A=Q-7 cosmology throughout this paper. 



2 RADIO DATA AT 610 MHZ 

We use the radio obse rvations of the X MM-LSS field de- 
scribed in detail in iTasse et al. I l|2007h . These observa- 
tions cover 13 degrees^ centered at a(J2000)=2''24"'00'' and 
5(J2000)=-4°09'47". This field was observed for a total of 
18 h with Giant Metrewave Radio telescope (GMRT) in Au- 
gust 2004 at a frequency of 610 MHz. The observations have 
a resolution of ~ 6.7" and a 5(j rms sensitivity r anging from 
1.5 m jy to 2.5 mjy across the observed field. iTasse et al.l 
l|2007h detect 767 radio sources in this field and obtain esti- 
mates or upper limits for the spectral indices of these sources 
by comparison with the NRAO VLA Sky Survey (NVSS). 



3 7f-BAND DATA 

We match the radio sources detected in the XMM-LSS 
field with _ftr-band sources detected in the 7th data re- 
lease (DR7) of t he ongoing UKIRT Inf rared Deep Sky Sur- 
vey (UKIDSS) (|Lawrence et al.ll2007^ . UKIDSS uses the 
Wide Fiel d Camera mounted on the UK Infrared Telescope 
(UKIRT) ICasah et al.l 120071 '). This survey, which began in 
May 2005, comprises 5 complementary sub-surveys of vary- 
ing depths and sky coverage. Our study uses data from 
the Deep extragalactic survey (DXS) which is a deep, wide 
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Figure 1. The location of the radio and infrared data used in 
this paper. The hght grey circles show the 610 MHz pointings 
of the GMRT. The thick black line indicates the UKIDSS DXS 
coverage of the XMM-LSS field in the 7th data release. 

survey over four observed fields. The survey is aiming to 
cover 35 degrees ^ to a planned 5a depth of K~20.8 (see e.g. 
iKim et al.l[2O10h . We use the online catalogue of sources in 
the XMM-LSS field, the sky coverage of this field in the 
7th data release is illustrated in figure [Uthe overlap in the 
DXS sky coverage and the radio observations is approxi- 
mately 2.75 degre es'^. This data is complete down to K=19.2 
(jKim et al.ll2010[ l and varies in depth across the observed 
field. Thus to ensure completeness we adopt a cut-off of 
K > 19.2 for the optical counterparts to our radio sources. 
K-band magnitudes in the U KIDSS survey are cali brated in 
the Vega system described bv iHewett et al.1 (|2006l l 



4 RADIO AND INFRARED MATCHING 

For the purposes of identifying near-infrared counterparts to 
the radio data we divide the radio data into two categories. 
The first category comprises sources whose radio and opti- 
cal emission are expected to be physically co-incident and 
includes all single component radio sources as well as par- 
tially resolved multiple component sources and double radio 
sources with clear radio cores. The second category contains 
sources where the optical emission is expected to be sep- 
arated from the radio emission, this includes double radio 
sources with no detected radio core, radio jets and other 
complex multiple component structures. 

To determine the optimal pairing radius for sources in 
the first category we consider the accuracy of the radio 
and K-haad source positions. The radio source positions 
in the ca talogue we use are accurate to approximately 2" 
ijTasse et al. 2007) and the UKIRT survey provide s accurate 
astrometry to within 0.1" l|Lawrence et al.ll2007l ). We thus 
adopt a pairing radius of 5", corresponding to roughly 2.5a, 
where a is the combined error in the positions of the radio 
and Tf-band sources. This large radius ensures relatively few 



true optical counterparts will be located outside our search 
area. However as the surface density of optical sources in- 
creases there is an increasing probability that an an optical 
counterpart may appear within our adopted search radius 
by chance. To determine the probability of a given opti- 
cal counte rpart being su c h a sp urious alignment we use the 
method of lDownes et al.l (|l986r ). This method calculates the 
probability P that an optical source could lie at the given 
distance from the radio source by correcting the raw pos- 
sionian probability for the number of ways that such an ap- 
parently significant association could take place. We reject 
all matches whose corrected probability value P> 0.05. As 
a further consistency check we plotted the contours of all 
the radio sources in DXS XMM-LSS field over the optical 
images and inspected the data by eye. 

Determining the optical counterparts of the more com- 
plex radio sources in our second category is much less suited 
to statistical methods. Thus we determined the best match 
for sources in this category by means of visual inspection. 
We find matches brighter than K=19.2 for 131 sources out of 
a total of 213 radio sources detected in the DXS XMM-LSS 
field. 



5 ESTIMATING REDSHIFTS 

We use a method of redshift estimation, developed by 
ICruz et al.l (|2007l ). which utilises the tight correlation be- 
tween Ti'-band magnitude and redshift for radio galaxies. 
This K-z relation has been investigat ed for a number of radi o 
galax y samples includi ng the 3CRR ijLillv fc Longaiijll9"83 ). 
6CE feales et al.lll997l) 6C* (iJarvis et al.ll2001al ) and 7CRS 
sample s (Wil lott et al.ll2003l ) and continues to at least z=3 
(jjarvis et al.ll2001ah . 

The method developed in ICruz et al.l (|2007l ) uses a 
model of the distribution of radio galaxies as a function of 
redshift together with the linear K — z relation to generate 
a Monte Carlo simulation of radio sources which populate a 
synthetic K — z diagram. The advantage of this method over 
simply applying a linear fit to the K — z relation is that it 
incorporates information about the scatter in this relation 
and is thus able to provide a measure of the uncertainty in 
the output redshift estimates. 

The distribution of radio galaxies and K — z relation 
used in our simulations were based on the combined dataset 
of the 3CRR,6CE,6C* and 7CRS samples. [Cmz et al.l(|2007l ) 
find the radio galaxies in this sample to be well fitted by a 
function of the form : 

n{z) = ylexp{-[ELoffli(logioz)T} (1) 

with polynomial coefficients : 

A=197.96; ao=-0.39; a2 = 1.00; 04=1.47; 

ai = 1.17; a3=1.83; 05=0.38. 

IWillott et al] (|2003l 'l quote the best fit to the K-z rela- 
tion in this combined sample as the second-order polynomial 

K{z) = 17.37 + 4.53 log^o z + 0.31(logio zf. (2) 

The Monte Carlo simulations assume Gaussian devia- 
tions about a mean T^-magnitude obtained from the K—z re- 
lation. A constant dispersion for these deviations is adopted 
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for the en tire redshift range. We use the dispersion ak ob- 
tained by ICruz et al.l (|2007l ) via a fit to the same dataset 
used to obtain the K — z relation. This dispersion was found 
to be afc=0.593. 

The highly populated synthetic K — z diagram gener- 
ated by the Monte Carlo simulations can be used to ob- 
tain a photometric redshift probability density function for 
a source of any given K-band magnitude. This is achieved 
by extracting the redshifts of all the simulated sources in a 
narrow interval about the K-band magnitude in question. 
A fit to the distribution of these extracted redshifts pro- 
vides the required probability density function. These dis- 
tributions are best-fit by a logjg-normal distribution with 
probability density function given by : 



p{z\K) = 



In (10)272 



:exp{- 



2a2 



} (3) 



where and o are the mean and standard deviation for 
normally distributed random variable logjg(2). The best- 
fitting estimate for z is defined as 

Zest = 10" (4) 

and the asymmetric 68% confidence interval is defined as : 

10"-" < Zest < 10" + " (5) 

The distribution of the estimated redshifts for all the ra- 
dio sources in our survey with if— band counterparts is 
shown in figure [21 This distribution is compared with the 
redshift d istribution predicted by t he SKADS simulations 
l|Wilman et al, 200^ : Iwilman etal] I2OIO. ) for a radio con- 
tinuum survey with flux limits comparable to the survey 
used in this analysis. 

This method of redshift estimation does not adjust for 
the possibility that optical host ga laxy luminosity and radio 
luminosity may be correlated (e.g. IWillott et aLll2003h . This 
would most likely occur as a result of a correlation between 
both thes e properties with the mass of the central black hole. 
lMcLure~et al. (2004) investigated the relationship between 
radio luminosity and black hole mass and concluded that 
these properties were only signiflcantly correlated for the 
higher luminosity FRIT and HERG sources, and the effect 
is relatively small. As we are predominantly interested in 
investigating the evolution of low luminosity sources in this 
analysis we are confident that this will not have a significant 
effect on our conclusions. 

Although using this method essentially precludes us 
from investigating any correlations between galaxy mass 
and radio emission, (e.g. IWillott et ai]|2003l : iMcLure et al.l 
I2OO4I '). it does allow us the possiblility of just using single 
band photometry to measure the evolution of radio sources, 
and is therefore much less telescope intensive than using 
multiband photometry. This is particularly true if one wants 
to constrain the brighter end of the luminosity function, and 
therefore driven to use larger survey areas. 



6 V/Vmax 

The evolution of the comoving space density of sources in a 
complete sample c an be assessed using th e non-parametri c 
V"/Vmax method (|Rowan-Robinsonl 1 19681 : ISchmidtl Il968h . 
Here V is the volume enclosed by the object and Vmax is 
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Figure 2. The distribution of the estimated redshifts of the 
matched radio sources in our sample. The top figure compares 
the best fit redshift estimates determined by our method to 
the theoretical redsh ift distribution in the SKADS simulations 
llWilman et al.l [2OO8I ) of a radio survey the same flux limits as 
that described in section [2] The bottom figure compares the best 
fit distribution to the distributions which would occur if all the 
sources had redshifts corresponding to the upper and lower 68% 
confidence intervals of their predicted redshift. 



the volume enclosed by the object if it were located at the 
redshift Zmax where it's measured flux drops below the limit 
of the survey. In order to determine Zmax it is necessary to 
take into account flux density limits of the radio survey and 
the limiting magnitude of the infrared survey. If Zmax de- 
termined by the radio data is larger than the redshift Zopt 
which would cause the source to drop out of the infrared 
survey then ^max = ^opt- We use our Monte Carlo simula- 
tions of the K-z diagram to predict Zopt- We obtain a value 
of 2opt = 2.5 as the mean z corresponding to our limiting 
K-band magnitude of 19.2. 

If the objects in our sample are uniformly disributed 
in space then V/Vi-na,x will be uniformly distributed between 
and 1 with a mean value of (V/Knax) = 0.5±(12Ai')- 2 , 
where N is the number of objects in the sample. A value 
of {V/Vimux} > 0.5 implies that comoving space density of 
sources increases at higher redshifts ie negative evolution 
with cosmic time whereas (V/Vmax) > 0.5 indicates a decline 
in space density with increasing redshift. 

The original {V/Vma.^} test was formulated for a sin- 
gle sample complete above a given radio flux density and 
optical magnitude. However the ra dio survey used i n this 
analysis is inhomogeneous in depth (jTasse et al.ll2007l l. thus 
it is mor e appropriate to use th e generalised 14/14 test de- 
vised bv lAvni fc Bahcalil |l980) to analyse the evolutionary 
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Figure 3. Average V/Vmax in seven different radio lum i nosity 
bins compared to the results obtained in IClewIev fc Jarvid l|20oi l 



properties of this sample. The generalised test treats the in- 
complete radio data as a single sample with a variable survey 
area. This variable survey area depends on the flux density 
of the individual radio source. The new test variables are 
the enclosed volume Vc and the available volume 14, the 
statistical properties of the Vc/Vk variable are identical to 
those of the V^/Knax v ariable in a single flux-limited sample 
ijAvni fc Bahcalll 1 19801 1. 

To determine the effective survey area as a function of 
radio flux density we perform a simple simulation. We in- 
ject a 1000 radio sources with peak flux densities ranging 
from 1-4 mjy into the original map of the radio data. These 
simulated sources are restricted to the 2.75 degrees^ cov- 
ered by the DXS survey. We determine the fraction of these 
simulated sources which are recovered by the source extrac- 
tion algorithm SAD in the Astronomical Image Processing 
Software (AIPS) as a function of their input flux densities. 



7 V/Vmax results 

In flgureOwe compare the mean K/Va statitistics obtaine d 
in our analysis with those found in lClewlev Sz Jarvij l|2004 ). 
We find that our results are in agreement with their anal- 
ysis which is based on galaxies selected at 325 Mhz in the 
Sloan Digital Sky Survey. O ur results also comp are very well 
with the F/Knax results of iTasse etU] (|2008li . Their anal- 
ysis is based on a combination of 325 MHz and 610 MHz 
radio sources in the XMM-LSS field, the data used in their 
analysis overlaps the DXS field studied in this paper so it 
is encouraging that we obtain a similar result. The strong 
evolution of the h igh luminosity radio sources detected by 
IClewlev fc JarvisI (J004) is also in evidence in our figure [3] 
although this result is at a lower statistical significance due 
to the small size of our sample, similarly their is a little ev- 
idence that the low power sources are evolving with cosmic 
time. The 610 MHz flux densities of our radio sources were 
converted to 1.4 GHz radio luminosities using the spectral 
index estimates obtained bv lTasse et all (|2007l l by compar- 
ision with the NVSS survey. 

As discussed in section [5l there is some uncertainty in 
the redshift estimates used in this analysis due to the in- 
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Table 1. Results of (Vc/Va) analysis is seven luminosity bins. 
The column entitled Icr (MC) quotes the uncertainty in our 
(y/Vniax) estimates due to the uncertainty in our redshift es- 
timates. The last column gives the mean z of sources in each bin. 




0.2 0.4 0.6 0.8 

Figure 4. The banded V/Vmax in two luminosity bins. The 
bin with Li.4GHz > lO'^^ W Hz~^sr~^ is represented by 
grey triangles and the black crosses represent the Z/i.4GHz < 
1025 w Hz-lsr-i bin. 



herent scatter in the K — z relation. To characterise the 
effect of this source of error on our 14/14 estimates we use 
Monte Carlo methods. In these Monte Carlo simulations we 
assign each source a redshift estimate determined by the red- 
shift probability density function derived from its' K-band 
magnitude. We use these simulated redshifts to rederive the 
(14/14) in the seven luminosity bins in figure [S] We repeat 
this process 1000 times to obtain the 68% confidence interval 
of the Vc/Va estimates. This confidence interval is indicated 
in table [1] . The size of this source of error is comparable to 
the statistical errors due to the sample size. 

To further our investigation of this evolutionary be- 
haviour we determine the banded 14/14 statistic in two 
luminosity bins. A high luminosity bin corresponding to 
sources with I/i.4ghz > 10^^ W Hz~^sr~^ and a low luminos- 
ity bin with Li.4GHz < lO^'^ W Hz~^sr~\ The banded test 
calculates the mean 14 /Vk for all sources with z > zo for a 
range of zq values. All evolution below zq is thus masked out 
of the analysis allowing us to determine whether the evolu- 
tion of these two populations of radio sources changes as a 
function of redshift. The results of the banded test are pre- 
sented in figure |4] The error in the banded test as a result of 
the uncertainty in our redshift estimates is also calculated 
using Monte Carlo methods and presented in table [7] 

The results in figure |4] indicate that there is little ev- 
idence of evolution in the low power radio sources out to 
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8 RADIO LUMINOSITY FUNCTION 

Another standard means of quantifying the cosmic evolution 
of radio sources is to compare the measured radio luminos- 
ity function (LF) of these sources at different cosmic epochs. 
Changes in the luminosity function $z (L) are usually mod- 
elled as being due to one of two simplified evolution scenar- 
ios, either pure density evolution or pure luminosity evolu- 
tion. In pure density evolution the luminosity distribution of 
the radio sources is considered fixed while the number den- 
sity of radio sources changes with redshift. In this scenario 
the measured LF at a given redshift $z(i) is related to the 
local luminosity function ^a{L) via the following equation : 



(6) 



In pure luminosity evolution the radio sources undergo a 
change in their luminosity with time, this evolution is pa- 
rameterized as: 



+ 



(7) 



We calculate the "I>z(-£/) of our sample in four r edshif t 
bins using the standard l/\max method of ISchmidtl (|l968l ). 
Thus the radio LF and its corresponding error is calculated 
for a particular luminosity bin as; 



j — \ 'max 



E 



(8) 



where the j index runs over all sources in the luminosity bin 
L-AI/ — >-L+AL. We adjust the Vmaa; values to take into ac- 
count the redshift limits of the four redshift bins for which 
the luminosity functions were calculated. The results of our 
analysis are presented in figure [5] and compared with the 
LF of radio sources in the VLA-COSMOS survey derived 



Table 2. Results of the banded (Ve/Va) test. The columns enti- 
tled la (MC) list the Monte Carlo estimates of the error in our 
(Vc/Va) estimates due to scatter in the K — z relation. 



a redshift of z ~ 0.8. Beyond this redshift there is a slight 
decline in the average Ve/Vk of this population however the 
errors derived from our Monte Carlo estimates increase sig- 
nificantly at these higher redshifts. The stronger evolution 
of the high luminosity sources can also be seen in this fig- 
ure. The mean Vc/Va for the low and high luminosity bins 
is calculated as 0.5 ± 0.03 and 0.57 ± 0.04 respectively. 
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Figure 5. The measured radio lumi nosity function in fou r red- 
shift bins compared to the results of ISmolcic et all ||2009| V The 
solid line repre sents the analytic fo rm of the local luminosity func- 
tion derived in lSadler et~all \200% . 



bv lSmolcic et al.l (|2009l ). Our luminosity function estimates 
appear to be in good agreement with the results of this pre- 
vio us study. 

ISmolcic et al.l (|2009[ ) characterise the evolution of the 
radio sources in the VLA-COSM OS survey by adopting the 
LF given by Sadler et all l|2002h as the best representative 
of their measured local luminosity function and using a least 
squares fit to derive the a o and a l parameters for pure den- 
sity and pure luminosity evolu tion of this LF. The analytic 
form of the radio LF given in ISadler et al.1 (|200 ^ was ob- 
tained from radio galaxies in the NVSS survey matched with 
optical counterparts in the 2d F Galaxy Red s hift Su rvey. 

The fitting procedure of ISmolcic et al.l (|2009l ) yeilded 
estimates of ao = 1-1 ± 0.1 and oil = 0.8 ± 0.1 where 
the errors were derived from the distribution of in the 
least squares fit. Employing this same method we obtain 
values for the density and luminosity evolution parameters 
of an ~ 0.6 ± 0.1 and aL = 0.8 ± 0.2 . These estima t es are 
fairly consistent with the estimates of lSmolcic et al.l (|2009l ) 
and imply an increase in the co-moving space density of low 
luminosity radio sources by a factor of ~ 1.5 in the interval 
between z—0 — 0.8. 

The con sistency of our results with those of 
ISmolcic et al.l (2009) provides encouranging evidence that 
the K — z relation can be used to investigate the radio lu- 
minosity function. This could provide a useful alternative 
technique for studying the LF in the absence of multi-band 
photometric redshift estimates. 
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Similarly ISadler et all (|2007l ) find evidence that low lumi- 
nosity sources experience mild evolution with an increase 
in their number density by a factor of ~2 at z=0.55. This 
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power AGN in the VLA-COSMOS survey out to z ~ 1.3. 
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As our analysis cate- 



gorises our radio sources solely on their luminosity we are 
unable to determine whether the separation in evolutionary 
behaviour is due to the FRI/FRII dichotomy or to differ- 
ences in black hole accretion modes in HERGS and LERGS. 

Our results demonstrate that using the /('—band (or 
similar wavelength), combined with radio surveys, is a viable 
route to investigating the evolution of the radio source popu- 
lation, at least up to z ~ 1.2. In the future, all-sky radio sur- 
veys such as those carried out wi th the Low- Frequency Array 
(LOFAR; iMorganti et alj|2010l ') and the Australian Square 
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(WISE), will enable us to pin down the evolution of the radio 
source population to a much higher degree of accuracy. 
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